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Abstract 

This is a written version of two lectures given at the First School on 
Field Theory and Gravitation in Vitoria (Brasil) April 15-19, 1997. 
We discuss the foundation of QCD as the theory of strong interactions 
and the perturbative analysis of e'''e~ annihilation to hadrons. Typical 
concepts of perturbative QCD studies, such as collinear singularities, 
jets, Sudakov form factors, are explained working out this case. 
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1 Introduction 



Quantum Chromodynamics (QCD), introduced by Gell-Mann and Frizsch in 
1972 JlSl, is the current theory of strong interactions. It is a renormahzable 



nonabehan gauge theory p9[ based on the group SU{3), containg quark and 
gluons as elementary fields, 

L= -}-F;^F^^^ + Y.^^{r^^D^-mf)qf (1) 
^ f 

where F^^ = df,A'^-d„A''^-gf''^''A''^Al is the field tensor, = d^. + igAp'' is 
the covariant derivative, t"" are the Gell-Mann matrices, f"''"^ are the structure 
constants of SU (3) and / denotes a quark fiavor. 

Strong interactions present many different phenomena and QCD at present 
cannot account for all of them. To give an example, QCD is irrelevant to the 
entire field of low energy nuclear physics. This situation is to be contrasted 
with QED, where most processes can be computed with high accuracy. 

In sec. 0we discuss the foundation of QCD, i.e. why we believe we have 
the correct theory (the lagrangian (P) even though many phenomena are 
not described or computed inside it. 

In sec. 1^ we discuss the fundamental property of QCD, asymptotic free- 
dom, according to which QCD approaches a free theory in the ultraviolet 
region. This property explains parton model assumptions and allows pertur- 
bative computations as corrections to the free quark behaviour. 

In sec. H a qualitative discussion is presented of quark confinement, ac- 
cording to which quarks do not exist as isolated particles, but only in bound 
states, the observed hadrons. This phenomenon cannot be described by 
definition in perturbative QCD, which assumes that asymptotic states, the 
hadrons, are infinitesimally close to the free ones, quarks and gluons. There 
are indications of confinement from lattice (non perturbative) QCD computa- 
tions, but a rigorous proof is missing. This section describes some qualitative 
ideas about confinement, to understand why perturbative QCD can be ap- 
phed despite confinement. We discuss also the production of hadronic jets, 
a typical phenomeon of high energy hadronic processes. The relation with 
confinement is explained, and also why we believe that a perturbative QCD 
computation can describe seminclusive jet properties. 

In sec. ^ we consider the production of hadrons in e^e^ annihilation 
at high energy as a typical application of perturbative QCD. This is a sim- 
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pier process than lepton-hadron or hadron-hadron scattering, so it is a good 
starting point for introducing the reader to perturbative QCD. The total 
hadronic cross section is perhaps the easiest thing to compute. Furthermore, 
the assumptions on hadronization to apply perturbative QCD are minimal: 
we require only that hadronization does not change the probability of the 
partonic process. 

Later we study the structure of the final hadronic state, which reveals 
many typical perturbative QCD phenomena, like coUinear singularities, jets, 
infrared effects to be factorized, etc. In the latter stronger assumption 

about hadronization has to be made: the momentum flow of the partonic 
event has to remain substantially the same after hadronization. 

There is also a section with the conclusions and an appendix describing 
the computation of e~^e~ annihilation into three jets. 

2 Foundation 

QCD is a non abelian gauge theory based on the group SU{3) of color. 
The main motivation for a non abelian gauge theory came from the parton 
model of deep inelastic scattering, according to which quarks behave like free 
particles in hard collisions. This model has been successively extended, with 
similar success, to describe also hard hadron-hadron collisions and hadronic 
production in e~^e~ annihilation at high energy (sec. ^). The hystory of deep 
inelastic scattering and the ideas of the parton model are sketched in sec. 

The choice of the gauge group and of the quark representation is moti- 
vated by various phenomena which are discussed in sec. |2.2| . 

We believe that all these different phenomena, collectively, give a good 
support to the structure of QCD as the theory of strong interactions. 

2.1 The Parton Model 

A series of experiments were performed at SLAC in the sixties to understand 
the structure of the proton and of the neutron, measuring the so called electric 
and magnetic form factors. Electrons with an energy up to 20 GeV were sent 
against a target of hydrogen or deuterium 

e + N ^e + X (2) 
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where X is any hadronic final state. 

The astonishing resuh was that a larger number than expected of large 
angle deflections of the electron were observed. Feynman gave a simple phe- 
nomenological explanation for this result: the nucleon has to be considered 
in deep inelastic collisions as a gas of non interacting pointlike particles, the 
partons; the electron simply suffers an elastic collision with a parton p 

e + p — > e + p (3) 

A pointlike cross section has not the form factor suppression of an extended 
object. We have therefore 'hard' interactions, with large angle deflections of 
the electron, as observed. To have a physical picture, consider the nucleon in 
a frame in which it has a relativistic velocity {'^Lorentz ^ 1), for example the 
rest frame of the colliding electron. Due to the Lorentz boost, the nucleon 
looks like a bunch of collinear partons, each carrying a fraction Xq of the 
nucleon momentum, such that 

= 1. (4) 

A hadron is represented by a function (partonic density) q{x) telling how 
many partons dn there are with momentum fraction x: 

dn{x X + dx) = q{x)dx (5) 

The observed hadronic cross section is the convolution of the parton density 
with the pointlike, partonic cross section. 

(Thadr = q*Crpart (6) 

Since the pointlike cross section assumes free partons, all the dynamical 
effects of strong interactions are contained in the specific form of q{x). We 
see here the idea of factorization at work: different dynamical processes are 
represented by separated factors in the cross section. This idea has been 
gathered and generalized by perturbative QCD, in many different forms: 
factorization of mass singularites. Operator Product Expansion (for a 
detailed discussion see for example Sudakov form factors (see sec. ^]3| 

and ref. [|T0[), etc.. 
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Parton model was also able to explain a phenomenological law, the so 
called Bjorken scaling, i.e the independence of the form factors on Q^, which 
coincide just with the parton densities. 

These partons were identified with the quarks introduced to classify the 
variety of hadrons discovered and their spectroscopy. Partons were assigned 
a spin s = 1/2, an electric charge e = 2/3, —1/3, a flavor, etc. That gave rise 
to a series of sum rules relating parton densities of different hadrons, which 
were experimentally satisfied. 

This identification between 'low energy' and 'high energy' entities came 
out to be non trivial. The valence quarks, the degrees of freedom of the quark 
model, were not sufficient to account for momentum conservation in e — iV 
collisions. The sum rule (^) seemed to be violated, the right hand side being 
substantially less than one. It was necessary to postulate the existence of 
'sea' quarks, in addition to the valence quarks, short living quark-antiquark 
pairs in the nucleon, partecipating to the hard processes. This idea was 
quite reasonable from the quantum field theory point of view, dealing with 
virtual particles form the very beginning, but was still incomplete. Half of 
the momentum of the nucleons was not carried neither by the valence quarks 
nor by the sea quarks because: 

^x + ^x ~ 0.5 (7) 

val sea 

It was necessary to admit the existence of neutral particles in the nucleon 
not interacting with the probe (the 7, W or Z), i.e. without electric or weak 
charge. These components were naturally identified with the gluons after the 
rise of QCD. 

2.2 The Color 

Nonabelian gauge theories are a good candidate to describe strong interac- 
tions because of asymptotic freedom (they are indeed a unique candidate in 
the framework of relativistic local field theories). The choice of the gauge 
group is motivated by a variety of phenomena, some of which are listed below 
(see also and the lectures of B. Mele at this School for a more recent 
compilation). These phenomena per se require only the color group as an 
(exact) global symmetry. Color enters as a multiplicity factor in the cross sec- 
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tion (counting of degrees of freedom) or as a quantum number to distinguish 
quarks. 



2.2.1 Spectroscopy 

Consider the A++ in the framework of a constituent quark model. It is a 
spin 3/2 baryon and is composed of three identical up quarks. Since it is the 
lowest lying state with these quantum numbers, it is natural to assume that 
the spatial wave function is the fundamental one. A 3-body potential which 
gives a reasonable description of the lowest excitations is the harmonic one: 

V{xi,X2,X3) = —^{{Xi-X2)^ + {XI-X3)'^ + {X2-X3)^) (8) 
D 

where Xi are the quark coordinates, ni is the constituent up quark mass and 

a; is a constant to be determined experimentally. 

The ground state of the model defined by the Hamiltonian 

H = -^t^ + Vi.u^,,.,) (9) 



IS a gaussian, 



ipo{xi, X2, xs) = — exp 



-({X,-X2f + {X,-X3f + {X2-Xsf)/{6R^ 



(10) 

where = l/{muj). ipQ is symmetric under exchange of any pair of coordi- 
nates and it has zero angular momentum (it is invariant under rotations): 

L = 0. (11) 

The spin of the baryon 

J ^ L + S (12) 
is enterely carried by the total spin of the quarks S: 

J = S. (13) 

The state with Sz — 3/2 consists of the up quarks all in the positive z 
direction 

|A++,5, = 3/2) = |«T,«T,wT) (14) 
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The spin wave function therefore is also symmetrical under exchange of the 
spins, so the complete wavefunction is symmetrical under exchange of the 
quarks. This is in contradiction with the spin-statistics theorem (the old 
Pauli principle of atomic physics) according to which identical spin 1/2 par- 
ticles must have an antisymmetric wavefunction. A natural solution of this 
problem is the introduction of an additional quantum number of the quarks, 
called the 'color'. There have to be at least three colors. The most eco- 
nomical solution is that quarks come in just three varieties of color, and the 
baryon wavefunction is antisimetrical under exchange of the color of any two 
quarks, so it is a color singlet. In color space the wavefunction is therefore 

I A++) = -^eabc\Ua,Ub,Uc), (15) 

where eabc is the totally antisymmetric tensor with €123 = 1. 

Similar problems exist also for the statistics of the A~ and Q' baryons, 
which are composed of three identical d (down) and s (strange) quarks re- 
spectively. Assuming SU{3) flavor symmetry, the same problem holds for 
many light baryons and is solved in the same way by the color. 

Mesons consist of quark antiquark pairs with the following color singlet 
wavefunction 

I M) = -^St I q") (16) 
where q and q' denote two arbitrary quark flavors. 



2.2.2 7r° 77 decay 

The 7r° has the following valence quark content: 

I = ^(1 I dd)) (17) 

We may take as an interpolating field for the 7r° (an interpolating field is an 
operator which excites the particle acting on the vacuum) the divergence of 
the isospin I — 1, 13 = axial current 

= -^{uJi^JbU-dj^j^d) (18) 
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Due to the triangle anomaly, there is a non zero coupling of the axial current 
with two electromagnetic currents, 

(0 I TA^{x)J,{y)Jp{z) I 0) 7^ (19) 

which produces the decay of the 7r° into two photons. The computation gives 

where /tt is a contant determined from leptonic charged tt decay (assuming 
isospin symmetry) and 

e = N{el ^ = 1 (21) 

with the standard charge assignements and N = 3. 

The decay amplitude is proportional to N because the pion couples to qq 
pairs of any color. 
The measured value is 

Tmis = 7.7±0.6e\/ (22) 

in good agrement with the = 3 rate. Note that this case is very sensitive 
to the number of colors, due to the quadratic dependence of the rate on N. 



2.2.3 Anomaly Cancellation In The Standard Model 

A pretty theoretical argument for the color (i.e. not based on any experi- 
ment) is the anomaly cancellation in the Standard Model. Anomaly has to 
cancel because a non vanishing anomaly induces terms of canonical dimen- 
sion in the lagrangian which are not originally present and are therefore not 
compatible with gauge invariance. The breaking of gauge invariance renders 
the theory not renormalizable, destroying the simplicity of the model (the 
simple ultraviolet structure). 

No triangle anomaly is generated if the sum of the charges for every 
generation is zero. Since 

e^ + ed = 1/3 (23) 

while 

e, + ez = -1, (24) 

a multiplicity for the quark states of 3 is required, which is provided by the 
color. 
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3 Asymptotic Freedom 



After the success of the parton model, the theoretical problem was that of 
recoinciling a good phenomenological model with field theory (so succesful 
in the area of electromagnetic and weak interactions). That appeared to be 
a hard task. Why was a quark behaving as a free particle when involved in 
an energetic collision? What about his interactions, which actually give rise 
to his binding in the nucleon? 

The resolution of this problem came with the renormalization of non- 
abelian gauge theories. In such theories the one-loop /3-function is negative. 
The /5-function is defined as 

Pi<^s) = /i^^ (25) 

and represents the variation of the renormalized coupling as we vary the 
renormalization point fi (see also sec. to understand its physical mean- 
ing). In this variation we have to keep the bare couplings (al?^ and the bare 
masses mo's) and the ultraviolet cutoff Aq fixed to ensure that we remain 
within the same physical theory (parametrized in many different ways ac- 
cording to the many different choices of the scale /z to which the parameters 
of the theory can be defined). Alternatively, we may define the /3-function as 

f^ia's) = A^^, (26) 

i.e. the variation of the bare coupling as we vary the ultraviolet cutoff Aq. 
The variation is done keeping the renormalized parameters constant, so as 
to remain in the same physical theory, i.e. to reproduce the same low-energy 
cross sections (with energy s <^ Aq). It can be shown that the two above 
definitions give the same function up to two loops [Q, beyond which details 
of the regularization (the precise way we cut the ultraviolet region) produce 
differences in the coefficients. The /^-function has a perturbative expansion 
which we write in the form 

(27) 
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We divided [3 by An because the one-loop beta function is already 0(a|) 
The one-loop value is 

= -11 + -nf = - 7.67 (/or Uf = 5) (28) 
3 

where iV = 3, is the number of colors (as proved in the previous section), and 
Uf is the number of 'active' flavors at the scale /x. With / < llA^/2 = 16.5 
flavors, i.e. / < 16, 

^0 <0 (29) 

The differential equation (^) can be integrated approximating the /3-function 
with its one-loop value. Denoting the initial integration condition with the 
scale /i^ and the final one with the scale Q^, we have 

' 1 - /3o/(4vr)a5(/x2) log(gV/i^) ^ ^ 

We see that, since Pq < 0, the coupling asiQ"^) decreases with the energy 
and approaches zero in the limit of a a very large energy (or momentum 
transfer): 

asiQ^) ^0 for Q2 ^ oo (31) 

The theory therefore approaches a free theory in the ultraviolet region (asymp- 
totic freedom ||25|), thus explaining the 'free' behaviour of the quark in the 
hard process. On the other side, the coupling is not small at low energies, 
explaining quark binding in the hadrons. We have therefore a very simple 
explanation of the duality of high-energy and low-energy hadronic phenom- 
ena with renormalization group ideas. We note that asymptotic freedom is 
a one-loop property of the /5-function because at large enough energies the 
coupling is small and higher orders are negligible. 

Instead of introducing a reference scale /z^ at which assign the value of 
the coupling, we can write the running coupling in terms of the so-called 
Aqcd parameter, the scale at which the one- loop coupling has a pole (it is 
an infrared pole because it occurs at low energies). With 

Al^^ = ^2gM(/3o<.,(M^)) (32) 
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the coupling is rewritten as 



1 



(33) 



-/?o/(4vr) log(QVA^eD) 

We may ask what is the physical meaning of Aqcd- It is a renormalization 
group invariant quantity (it does not change if we vary /z^ in its defining 
equation) and denotes at which energy renormalization group improved per- 
turbation theory ceases to be valid, pointing therefore to generic nonpertur- 
bative phenomena This scale has to be determined experimentally and 
turns out to be ~ 300 MeV, i.e. of the same order of the hadron masses 
{Mp ~ 770 MeV). 

We note that the coupling goes to zero only logarithmically with the 
energy, i.e. very slowly. As long as <^ 1, perturbative corrections to the 
free theory behaviour can be computed and are expected to be sizable. This is 
the technical justification for perturbative QCD, whose physical justification 
will be discussed later. 

Let us quote also the values of two and three loop coefficients of the 
/^-function in the MS scheme [0: 

Pi = 



38 

-102 + —Uf 
3 ^ 



P2 



1 r 
2 



2857 



5033 
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-Uf 



102 + 12.66 Uf 
325 



38.66 



27 ^ 



-1428.5 + 279.6nj - 6.0n- 



/ 



- 180.5 



(/or Uf = 5) (34) 

We note that /5i < if nj < 8, while < if < 5. The (approximate) 
two-loop solution for the running coupling is given by: 



-/3o/(47r)log(gVA2 



QCDJ 



1 + 



iog(iog(gVA2 



log(QVA2 



In QED the /3-function is given to one loop by 



Pqed 



3^ 



> 



(35) 



(36) 



(times Hf if f flavors are considered). Integrating the related differential 
equation exactly as in the case of QCD, we derive: 
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where A/,, called the Landau pole, is now an ultraviolet scale ||T9[. The 
coupling grows with the energy as a consequence of the positive /5-function. 
The QED behaviour of the effective coupling is therefore opposite to that of 
QCD (eg. (p3D). A physical discussion of these different behaviours is given 
in the next section. 

We conclude this section observing that in four space-time dimensions 
asymptotic freedom is an exclusive property of non abelian gauge theories 

which therefore emerge as a natural candidate for the strong interactions. 

3.1 Antiscreening 

In QED the behaviour of the effective charge is controlled by the photon 
polarization diagram, 

7* ^ (e+e-)* ^ 7* (38) 

(ultraviolet divergences in the electron field renormalization and in the vertex 
correction cancel because of the Ward identity). 

The physical explanation is the following. Consider a charged particle 
placed in the vacuum, with bare (i.e. primordial) charge cq. Quantum field 
fluctuations induce short lived e"*"e~ pairs. These pairs constitute electric 
dipoles which align with the test charge so as to minimize the energy. The 
bare charge is therefore surrounded by a shell of opposite charges (the outer 
shell is at infinity) with a consequent screening effect. The vacuum behaves 
as an ordinary medium with material electric dipoles, with dielectric constant 

e > 1 (QED) (39) 

If we measure the charge of the test particle, we find that it is less than the 
bare charge: 

e = ^ (40) 

If we come closer to the test particle, say at distance r, we see an increasing 
charge, because the outer part of the shell, at distance r' > r, does not 
contribute to the screening (Gauss theorem). The effective charge (i.e. the 
observable one) is therefore a function of the distance 

e = e(r). (41) 
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It increases as we lower the distance, or equivalently, as we increase the 
momentum transfer ~ 

In QCD the situation is reversed: due to gluon-gluon couplings (three 
and four gluon vertices), the test charge is surrounded by gluons of the same 
charge. The vucuum behaves as an hypotetical medium with dielectric con- 
stant 

e < 1 {QCD) (42) 

so we have an antiscreening effect. As we come closer to the test particle, 
the effective charge is reduced (see for a more detailed explanation). 

Asymptotic freedom means that the measured charge approaches zero at an 
infinitesimal distance. 



4 Confinement 

Isolated quarks have never been observed, as particles with fractional electric 
charge, for example. Hence the hypotesis that they are 'confined' in hadrons, 
i.e. that they do not appear as asymptotic states in any physical process. 
The current explanation is that the force between two quarks does not vanish 
as their separation r increases, so that the potential energy of the system V 
diverges 

V{r) oo as r ^ oo. (43) 

In any realizable collision, the available energy is finite, making therefore 
impossible to produce isolated quarks. 

A potential which is belived to represent the QCD interaction between 
a quark and an antiquark in a color singlet state, is the so called 'funnel' 
potential: 

V{r) = --^ + kr (44) 
where is a constant, called the string tension. 

This potential has been extensively used in non relativistic models of cc 
and bb bound states, with good results, and is also compatible with lattice 
(non perturbative) QCD computations. We see that at small distances, r -C 
1/k, the potential is coulombic, as it is in perturbative QCD, while at large 
distances there is a linear grow of the energy with the separation. 



The qualitative explanation of the confining potential (|4^) is the follow- 



ing. When the separation of the quark-antiquark pair is small compared to 
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the confinement radius, the gluon field has the same form of a dipole field 
in classical electrodynamics. Chromoelectric field lines spread out in all the 
space (and the potential decays as 1/r like in QED). When their separa- 
tion becomes of the order of the confinement radius, there is a reciprocal 
attraction of the field lines, which makes them to collapse in a tube along 
the qq line. Since the electric fiux is the same through any closed surphace 
containing the quark or the antiquark, and the electric field is non zero only 
inside the tube, we say that there is the formation of a fiux tube (or string). 
A further increase of the separation does not change the form of the gluon 
field anymore but only makes the fiux tube longer. Since the gluon field 
is now essentially concentrated in a tube, the energy contained in the tube 
is proportional to its lenght. In other words, there is a constant energy of 
the field per unit lenght, which implies a linear term in the potential. In 
the next section we discuss a model of confinement which induces the above 
postulated attraction of field lines. 

4.1 Dual Meissner effect 

A model which induces fiux tubes (and therefore confinement via a linear 
potential) is the so called dual Meissner effect p. 

Consider a superconductor in an external magnetic field B above the 
transition temperature Tc, for definiteness a cilinder with its axis parallel 
to B. As we lower the temperature to T < Tc, there is the transition 
from the normal to the superconducting phase (drop to zero of the electrical 
resistence). At the same time, the magnetic field lines are expelled outside 
the sample. The magnetic field inside the cilinder is zero: 

Binside = 0, (45) 

i.e. the superconductor shows a perfect diamagnetism 

Xm = 0. (46) 

The explanation is that with the transition electrons pair in particles with 
charge 2e (Cooper pairs). The latter induce a current Jpairs along the sur- 
phace of the cilinder generating a magnetic field inside the sample opposite 
to the external one: 
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Now imagine of placing a pair of opposite charge magnetic monopoles inside 
the superconductor. Since the superconductor tends to expel the magnetic 
field, the field lines of the magnetic dipole are shrunk in a small tube con- 
necting the monopoles. An electric current of Cooper pairs develops along 
the surphace of the flux tube, effectively squeezing the magnetic dipole field. 

As well know, clcctromagnctism is invariant under the simultaneous ex- 
change of magnetic fields with electric ones and of magnetic charges with 
electric ones (duality symmetry). 

Imagine the dual of an ordinary superconductor, i.e. a system expelling 
electric fields. If we place a pair of opposite electric charges inside the sys- 
tem, the lines of the electric dipole field will be shrunk into an electric flux 
tube. The dual superconductor contains magnetic monopoles, while the or- 
dinary superconductor contains 'Cooper pairs, i.e. 'electric monopoles'. Let 
us describe the phenomenon a little bit more in detail. Soon after the in- 
troduction of the electric charges, the system reacts producing a current of 
magnetic monopoles Jm along the surphace of a tube connecting the test 
charges. Jm produces, by duality, a secondary electric field Em cancelling 
the external field E outside the tube 

Em = -E (outside) (48) 

and reinforcing the external field inside the tube. The net effect is the squeez- 
ing of the dipole electric field of the pair into a tube. 

Confinement has not been rigorously proved in QCD up to now. It is 
usually assumed that confinement means that hadrons consist of color singlets 
only. 

4.2 Jets 

The process of jet formation in high energy collisions is qualitatively related 
to confinement. As for confinement, there is no rigorous proof of this phe- 
nomenon and of the ideas explaining it inside QCD. 

Let us consider the simplest case of jet production, e~^e~ collisions, to 
be discussed in detail in the next section. The electron and the positron 
annihilate into a photon, which decays into a quark-antiquark pair: 

e'^ + 6" 7 q + q. (49) 
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In the center of mass frame, the quark and the antiquark have the same 
energy and opposite spatial momenta (back to back). They are originally 
created in the same spatial point (the fireball) and fly far apart with light 
velocity. When their separation comes close to the confinement radius, 

Tconf ~ 1 /m = 10"^^ cm, (50) 

the dipole color field of the quark pair shrinks to a tube (called flux tube 
of string). As they separate from each other, the flux tube gets longer and 
an increasing fraction of the kinetic energy of the pair is converted into field 
energy. When the energy contained in the flux tube exceeds the mass of a 
light qq pair 

Etube > 2m, (51) 

a novel quark-antiquark pair is created which screens the color of the original 
pair. The string breaks into two shorter strings, each having at his ends an 
original quark and a novel quark in a color singlet state. 

If the relative momentum of the quark-antiquark pair connected to the 
same string is large enough, the string may break again. We have therefore 
the production of a parton cascade, i.e. of a large number of quarks and 
gluons with progressively lower energy. 

In high energy processes {E 3> Mp, say, where Mp is the proton mass), 
the flux tube has a large longitudinal momentum {Pl ~ E), while it has a 
limited tranverse momentum (P^ ~ Mp). Consequently, the partons pro- 
duced have large longitudinal momenta and large longitudinal momentum 
differences, but limited transverse momenta. 

After the parton cascade production, we have the recombination of quark 
and gluons in color singlet states, the observed hadrons. This transition 
does not induce large momentum transfer between the partons (it is a soft 
process) . This hypotesis is sometimes called local parton-hadron duality and 
states that hadrons are produced by partons which are close in phase space 
[ P^ . This mechanism forbids, for example, that a quark with momentum 
p coming from one jet combines with a quark with momentum p' coming 
from another jet {p ■ p' ^ Ag,^^), destroying the shape of the perturbative 
momentum distribution. The transition from a partonic jet to a hadronic one 
therefore does not wash out all the parton information and the longitudinal 
and transverse momentum distributions are substantially unchanged. The 
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final particles, the liadrons, are collimated into two small angular regions 
(jets) opposite to each other see (fig.l). 

Furthermore, the angular distribution of the jets relative to the e^e~ 
flight direction is expected to coincide with that one of the original qq pair. 
As we will see, this prediction is confirmed by data, giving good support to 
the above qualitative ideas. 

If an energetic gluon is produced at large angle with respect to the original 
pair (a hard gluon), a new string is formed with a consequent third jet (see 



From these considerations, it emerges that jets are a universal occurence 
of high energy hadronic reactions. 

5 Hadron production in e+e^ collisions 

Another evidence for the color comes from the hadronic production in e~^e~ 
machines (in this section we rest on ref. |]12|)- Consider the production of a 
fermion pair In lowest order, the electron and the positron annihilate into a 
7 or a Z^, which decay into a fermion- ant if ermion pair (/ 7^ e), (see fig. 3) 



fig.2). 



+ e 




/ + / 



(52) 



The differential cross section is given by: 



dcosO 



da 




where 9 is the angle between the initial and the final pair. 




s{s - Ml) 



(54) 



k 



1.40 



(55) 



Ana 
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{k is the adimensional ratio of the relevant weak couphng to the electromag- 
netic one and its numerical value is for a{Mz) = 1/128 ^^). Gp is the Fermi 
constant, a is the fine structure constant, Mz and Tz are the mass and the 
width of the Z boson respectively. Vf and af are the vector and axial vector 
couplings of the Z to the fermions and are given by: 

Vf = L^f — 2efsm'^ 6w 

ttf = hf (56) 

where 9w is the Weinberg angle, and J3 is the z-component of the weak 
isospin: J3 = 1/2 for neutrinos and u-type quarks (e = 2/3), and /a = —1/2 
for charged leptons and (i-type quarks (e = —1/3). 

The function Xi{s) represents the interference of the 7 amplitude with 
the Z one, and the function X2(s) represents the Z amplitude squared. 

At small energies compared to the Z mass 

s < M| (57) 
the functions Xi ^"^^ X2 are suppressed by the Z mass: 

xi ~ -k^ 1, X2 ~ X? < 1 (58) 

Z 

We can therefore neglect both the direct Z effect and the interference one 
(Xi = X2 = 0), ending up with the pure QED cross section 

_ = _e}-(l + cos^9) (59) 

The angular dependence is in good agreement with the observed 2-jet cross 
section and gives support both to the spin=l/2 assignement of the quarks 
and to the qualitative ideas about jet formation discussed previously (see 
fig-4). 

LEPl and SLAC machines have operated on the Z peak, i.e. with a 
center of mass energy E = Mz- In this case the resonant contribution of 
the Z dominates over the 7 and the interference ones. Setting therefore 
6/ = 0, Xi = and X2 = ^^^|/r| we have for the differential cross section 



1 + cos o){a^ + v^)[aj: + Vj) + 8 cosuaeVeajVf 



dcos9 2T% 



(60) 
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The angular distribution contains an additional cos 6 term, produced by 
the interference of the vector and the axial current. This term induces a 
difference in the number of events in the hemispheres on the electron and 
the positron side {9 : (— 7r/2, 7r/2), (7r/2,3/27r) respectively). This effect is 
called forward-backward asymmetry and is generated by the C /P violation 
of weak interactions. 



5.1 The Ratio R 

Integrating over the polar angle, we derive the total QED cross section: 



a — 

Let us define the ratio 



a^e)- (61) 

6 ' S 



„ a(e^e —>■ hadrons) 

R = \ + _ 62 

This way we get rid of uninteresting factors and concentrate on strong inter- 
actions only. On the experimental side, some systematic effects, related to 
the luminosity determination for example, cancel in taking the ratio. 

We assume now that the total hadronic cross section is equal to the total 
partonic cross section, i.e.: 

a{e^e^ — > hadrons) = a{e^e~ — > partons) 

= a{e~^e~ qq) + a{e^e' qqg) + ... (63) 

There is a physical justification for that. The time for the partonic process 
to occur (virtual photon decay) is 

- ^ (64) 
where E is the center of mass energy, while the time for hadronization is 

th ~ -r^ (65) 

J^QCD 

(or equivalently, the lifetime of a resonance r ~ 10~^^sec). In high energy 
collisions, 

E > Kqcd, (66) 
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hadronization occurs therefore much later than the partonic process, 



tp < th, (67) 

so its cross section is substantially unaffected. The above qualitative expla- 
nation is not rigorous, but, as we will see, its conclusions turn out to be rigth. 
We have therefore in lowest order: 

R = N^Y.^) (68) 

/ 

where the sum extends to all the quark flavors which are kinematically al- 
lowed. The ratio R has therefore the value: 

R = N{el + el + el) = 3(4/9 + 1/9 + 1/9) = 2 (69) 

below the charm- ant icharm production threeshold, i.e. below a center of 
mass energy of ~ 3 GeV, a value of 

R = N (el + el + el + el) = 3 (4/9 + 1/9 + 1/9 + 4/9) = 3.33 (70) 

between the charm and the beauty production threshold, i.e. between ~ 4 
GeV and ~ 10 GeV, and 

R = N{el + el + el + el + el) = 3(4/9 + 1/9 + 1/9 + 4/9 + 1/9) = 3.67 (71) 

above the beauty threeshold and well below the Z peak (see fig. 5). The 
above estimate is basically a counting of the degrees of freedom involved. It 
is remarkable that it depends only on the electric charges and the colour of 
the quarks. The above prediction for R is valid only far away from cc and 
bb resonances, where hadronization effects cannot be considered negligible 
because they affect substantially the recombination of partons. We see that 
R in the resonance regions has peaks and valleys. In these regions we may 
interpret the QCD estimate ( pE| ) as an average over an interval of energy AE 
of order 1 GeV 

1 rE+AE/2 

R ^ IR) = R(E')dE' (72) 

^ ' AEJe-AE/2 ^ ^ ^ ^ 

In the low energy range (a few GeV), the parton estimate is not very accurate 
because hadron production is still affected by light resonances. Furthermore, 
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perturbative corrections (to be discussed later) are large, and strongly de- 
pendent of unknown higher orders. 

At E = 34 GeV the experimental value of R is about 3.9, to be compared 
to 3.67 (see fig.5). The mismatch is only reduced by the Z contribution (it is 
a 0(1%) effect, as can be evaluated inserting the low energy approximation 
for Xi ill eg. (|53D ) . As we will see, radiative corrections enhance the value of 
the lowest order, bringing the theoretical value in the experimental band. 

Let us consider now the resonance value {E = Mz) of R. Integrating over 
9 we have for the total cross section: 

<^ = ^-^{al + vl){a) + v)) (73) 

We note however that real photon emission from the initial e^e~ state (called 
bremsstrahlung) smears the peak and shifts its position (see fig. 6). So, the 
previous formula has to be corrected for an accurate quantitative analysis 



(see the end of sec. ^]3| for a discussion of these effects). The ratio R is given 
by: 



E f-1 flf + v'] 

R = 3 ^ , f = 20.09 74 

where we have used eqs.(|56|) and siv? Ow = 0.2315 We note that the 

QED effects mentioned previously (bremsstrahlung) are the same for the 
muonic and hadronic channels because they can be factorized as part of 
the initial e^e~ annihilation process. So they cancel in taking the ratio of 
the cross sections. The prediction for the tree level value of R on the peak 
therefore turns out to be accurate in this respect. Note that R is much larger 
than in the QED case because quarks have comparable weak charges to that 
of the /i. 

The experimental value is [^: 

Rexp = 20.80 ±0.035 (75) 

Also in this case, the theoretical value is ~ 4% smaller than the experimental 
one, because of corrections to be discussed in the next section. 

5.1.1 Radiative corrections 

We have seen in the previous section that the measures of R are so accurate 
that a comparision with the theory requires the inclusion of perturbative 
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corrections. Let us consider therefore radiative corrections to R. It is a com- 
pletely inclusive quantity and, neglecting quark masses, it is characterized 
by a single energy scale. 



(76) 



(quark masses can be sent to zero because in this limit no singularities are 
generated). The order as corrections are ultraviolet finite because the elec- 
tromagnetic current which creates the qq pair is conserved in QCD (see the 
appendix). In other words, the UV singularity of the vertex correction is can- 
celled by the UV singularity of the external quark legs renormalization. In 
the intermediate stage of the calculation, infrared divergences appear, which 
cancel between real and virtual diagrams in the completely inclusive process. 
In this case they offer therefore only a technical problem. Let us see this 
explicitely in dimensional regularization. Infrared and coUinear singularities 
(see later) are regulated increasing the space-time dimension to n > 4 (see 
for example |2^). The real corrections, with a gluon in the final state (see 
fig. 7), give: 



mg) 



CTq- 



27T 

27r 



H{e) 



dxidx2 

2 3 
^ + - ^ 



xf + xl — e(2 — xi — X2) 



(1 
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T 



-xiY+'ii 

+ 0(e) 



(77) 



where Cf = E tata = {N^ - l)/2N for an SU{N) gauge theory (Cp = 4/3 in 
QCD) is the Casimir operator in the fundamental representation, n = 4 — 2e 
and 



Hie) = 



3(1 



(3 - 2e)r(2 - 2e) 

The virtual corrections, i.e. the diagrams with a reabsorbed gluon, give: 



(78) 



a(e e 



Zn 



3 

e 



8 + 0(e) 



(79) 



Soft singularites appear as poles in e. The double pole 1/e^ is originated 
by the product of the coUinear with the infared singularity (see later), while 
the simple poles 1/e come from a coUinear or an infrared singularity. Both 
double and simple poles cancel in the sum, giving a finite 0(0:5) correction 
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to the tree level value of R: 



R = Rtree(l + —) (80) 

In the computation to order a|, we hit (no more cancelling) ultraviolet di- 
vergences, of the form 

«|log^, (81) 

where Aq is an ultraviolet cutoff, i.e. an ultraviolet regulator, and as = a^s^ 
has to be interpreted as the QCD bare coupling. To be more explicit, the 
correction factor to i? up to 0(a;|) has the form: 



K = l + C as-C ^al log ^ + C'ai + ... 

47r ^ s ^ 



2 l^^'^O , r^l2 



= l + C as{l-^as\og^) +C'al + ... (82) 

(as we have seen C = 1/vr). We can compare this result with the one-loop 
running coupling of QCD (eq.(|5^)), in which we set /i^ = Aq and we interpret 
^^(Aq) as the bare coupling (we also set = s): 

I . «s(Ag) 
as[s) - 



l + /?o/(4vr)as(Ag) logAg/s 
= <^s[l-j^aslog-+^-j a, log - + ...j (83) 

We may write, up to order q;|: 

K = 1 + Cas{s) + C'as{sf + ... (84) 

We see that the logarithmic term in K is simply a renormalization of the 
bare coupling. The theory is trying to say that we are expanding around a 
wrong scale, the ultraviolet cut-off Aq, very far from the physical scale s. The 
unfactorized cross section (^) and the factorized one (0) coincide to order 
a| included. The terms added to generate the factorized form are indeed of 
higher order: a| or more. 
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We see however that the above factorization is physically reasonable but 
not unique. We may equally well write: 



K = 1 + Cas-C^al log ^ + 



An 



a's + ... (85) 



with a consequent factorized formula of the form: 



K = 1 + Casis/A) + 



An 



asis/AY + ... (86) 



We see that the running coupling is evaluated at a different scale, s/4, instead 
of s, and there is a compensating term of order a|. So, there is no way in 
perturbation theory to fix the scale in the running coupling to a unique 
value. We can take whatever value for the scale, and compensate shifting the 
finite 0(a|) term. Of course, a wird scale setting of lO'^^s is unreasonable, 
because it introduces a large logarithmic term 4 log 10 ~ 10 in the finite a'g 
term deteriorating the convergence of the expansion. 

The two factorized formulae ( plj) and (p6| ) agree formally up to order a| 
but produce slightly different numerical values. They differ in the way they 
capture parts of the higher order terms. This phenomenon occurs because of 
the different treatment of the coupling as and of K. The running coupling 
resums classes of terms of any order in as, as it is clear from eg . (|8^) . The 
group structure of scale tranformations (the renormalization group) allows 
this resummation, while there is not any known way to resum terms in K to 
any order in as'- K is expanded up to a fixed order in as, i.e. the series for 
X is a truncated one. When we change the scale of the coupling, an all-order 
change, K cannot compensate for the tail of the series. The independence 
on the (arbitrary) scale choice is therefore not rigorously true in perturba- 
tion theory because it is violated by higher orders. This ambiguity never 
disappear: it can only be shifted to higher orders as we compute more and 
more terms for K. We have the formal independence from the scale choice 
of the coupling only in the limit in which the infinite series for K has been 
computed. 

Let us rephrase these facts in a more common language. The UV diver- 
gence (|8lD is converted with traditional renormalization in a similar logarith- 
mic term: 

a|log4, (87) 
A* 
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where /x is a renormalization point. In dimensional regularization (the only 
regularization which is effectively used in high order computations), after the 
subtraction of the 1/e poles, we are left with logs of the above form, where n 
is the mass unit introduced to make the coupling dimensionless in dimension 
n 7^ 4. 

R is therefore written as a function of the renormalized coupling at a 
fixed scale /x^, asifJ^"^), and of the above logs: 

R = Riasi^i)^ogs/^I^) (88) 

since i? is a physical quantity, it cannot depend on the choice of the renor- 
malization scale /z^, which is an arbitrary scale introduced for convenience 
at an intermediate stage (at the very end we compare R with another pro- 
cess characterized by another physical scale s'). These considerations are 
formalized with the renormalization group equation 

,^jLn . 0. (89) 

where as usual, one has to vary fi^ remaining in the same physical theory, 
i.e. keeping fixed the bare parameters of the theory, or the renormalized 
parameters at some other fixed scale /Iq. Explicitely eq.(PD|) reads: 

fi 7— + (3[as)- 



djj'^ das 



Rias,\ogs/fi^) = (90) 



The concrete prediction we are able to give is slightly dependent on the renor- 
malization point, an arbitrary scale which has neither a physical meaning nor 
a definite value. We need some physical intuition to pick up a good value of 
/i^ which minimize the (unknown) higher orders: a natural choice seems to 
be fj,"^ = s, for which i? is a (truncated) series in as evaluated at the center 
of mass energy. 

R = R{a{s)) (91) 

We have therefore: 

R = K*Rtree (92) 

where the correction factor (the so called i^-factor) has the following series 
expansion 

K=l + tcJ^)" (93) 
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In the MS scheme (an off-shell renormalization scheme, specific of dimen- 
sional regularization) , three correction terms are known 



Cl = 1 

c. - f -UC(3)+(fc(3)-H))„, 

= 1.986 - 0.115 nf = 1.411 

87029 1103^,^, 275^,^, tt^ /55 5^,^A 

C3 = C(3) + C(5) /9n+ C(3) h 

^ 288 4 ^ 6 ^ 48^° \72 3^^ 7 ' 

/7847 262,, , 25,, A /151 19,, A . 
= -6.637 - 1.2n/ - 0.005nJ - 1.240?; 

= -12.80 (94) 

where the last values of the coefficients are for nj = 5. ({s) is the Rieman 
C-function, defined by the series 

oo -1 

C(-) = E-. (95) 

n=l 



and 

7^ = = ^ = 0.0303 (96) 



3 E/e} 33 



for the electromagnetic diagram (the last value is for rij = 5) and 
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for the weak diagram. The different dependence on the electric/weak charges 
of the 7] terms is generated by diagrams with two separated fermionic traces 
(like for example in the process e^e~ 3g). 

At = 34 GeV, taking a,(34 GeV) = 0.146 ± 0.030, we have K ^ 1.05, 
i.e. the required 5% increase of the tree level value for R. At LEPl energies, 
y/s = Mz, taking as{Mz) = 0.12 [g, we have K ~ 1.039, i.e. a 4% increase 
of the tree level value. The experimental measure of R is therefore a clean 
verification of QCD. Note the large (negative) coefficient of the a| term. 
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which makes the perturbative series not so well convergent. For a<j = 1, for 
example, the q;| term is almost three times bigger than the a;| term. 

The dependence of the coefficients Cj on the scale, shown explicitely in 
lowest order at the beginning of this section, can be derived in general im- 
posing the Callan Symanzik equation (|90D order by order in as- 



5.2 Jets in e+e collisions 

Perturbative QCD predicts also some properties of the structure of the final 
state. Of course, the cross section for a single exclusive channel, like for 
example e^e~ — >■ PP, involves detailed hadronization processes and is there- 
fore outside the reach of perturbative QCD. We can consider seminclusive 
quantities. Let us consider the emission of a gluon, i.e. the process 

e+ + e" q + q + g (98) 

There are two amplitudes at order as, for the emission of the gluon from the 
quark leg Ma or the antiquark one Mb (fig. 7). The amplitudes involve only 
the QED-like vertex 

^mg = igi^ita, (99) 

i.e., the non abelian nature of QCD (three and four gluon couplings) does not 
fully manifest itself. The only effect of color is in the factor Cp multiplying 
as- This implies that the QCD cross section is identical to the abelian (QED) 
one with the replacement 

Cpas a (100) 
The non-abelian nature of QCD enters instead explicitely in the process 

e+e~ qqgg, (101) 

which is the dominant contribution to the 4-jet cross section and starts at 
order a|. Let's go back to the 3-jet case. We denote with Xi the energy 
fractions of the quark, the antiquark and the gluon: 

^1 = ^, X2 = ^, X3 = ^, (102) 
J^b -c/fe 
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where Eb is the beam energy {^/s = 2Eb). The conservation of energy and 
momentum gives: 



Xi + X2 + X3 = 2 

1 ~Xi = -XjXk{l - COS 6 jk), i ^ j k, i ^k. (103) 

The kinematical bounds are therefore 

< Xi < 1 (104) 

and 

Xi + Xj >1, i ^ j (105) 

In words, a parton has the beam energy as the maximum energy when the 
other two partons are colhnear, and a parton pair has the beam energy as 
the minimum energy and the total energy as the maximum energy. 
Since the final state is planar, we have also: 

^12 + ^23 + ^13 = 27r (106) 

The differential cross section for having a quark with an energy fraction xi 
and an antiquark with an energy fraction X2 in the final state is 

Ida^ = C,^--d±A-- (107) 
a dxidx2 27r (1 — a;i)(l — 0:2) 

This cross section is derived in the appendix; it coincides with the integrand 
in eq.(|77D in the limit e ^ 0. It is symmetric under exchange of Xi with X2, 
as it should be for the C/P symmetry of QCD (to order as, we divide da 
simply by ctq). 

If the experiment does not distinguish between the jets formed by the 
quark, the antiquark and the gluon, it is natural to symmetrize the above 
cross section: we simply ask what is the cross section for having a jet with 
energy fraction Xi, a second jet with energy fraction X2 and a third jet with 
energy fraction X3. The cross section is therefore: 



1 dasjets _ f xj + xl _^ xj + x"^ xl + x^ 
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a dx\dx2 Gtt Y (1 — xi)(l — X2) (1 — xi)(l — X3) (1 — X2)(l — X3) / 

(108) 
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Since the jets are assumed to be undistinguishble, The energy fractions Xi 
can be permuted so that 

Xi < X2 < X3 (109) 

We therefore say that jet 1 is the less energetic one, jet 2 the intermediate 
one, and jet 3 the most energetic one. Let's go back to the simpler formula 
( |107|) . The integration region for the total cross section is 



< Xi,X2 < I 

1 < X1 + X2 (110) 

There are singularities when the quark energy fractions reach their maximum 
allowed values (end point singularites) : 

Xi^l- i = l,2 (111) 

The strongest singularity is a double pole and occurs when the energy frac- 
tions reach simultaneously their end point values: 

Xi 1^ and X2 1^ (112) 

This happens when theee body kinematics 'resembles' the two body one (the 
latter is a double spike in the energy fractions ~ 5(1 — xi)5(l — X2))- 

As we have seen, the double pole ( p.l2| ) is converted with the integration 
in dimensional regularization in a double pole 1/e^ (eg. (|77|)). As we will see, 
the configuration ( |112| ) generates a double logarithmic contribution in the 



3-jet cross section, i.e. a term of the form log^ y. 
There is instead a simple pole singularity when only a single energy fraction 
Xi reaches the singularity, 

Xi 1", Xj < 1, j j^i (113) 

(by Xj < 1 we mean less than one by a finite amount, Xj = 0.5 for example). 
This is explicitely seen decomposing the kernel as 

xj + xl _ 2 1 + a;i 1 + X2 ^^^^^ 



(l-xi)(l-X2) (l-xi)(l-X2) 1 - a;2 1 - xi 

The configuration ( |113| ) is related to simple 1/e poles in dimensional regu- 
larization and to single logarithmic terms in the 3-jet cross section, i.e. to 
terms of the form as logy. 
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Both the leading and the subleading singularities produce divergences in 
the total cross section, so that 

crtot(e+e" qqg) = oo, (115) 

which is an unphysical result. Let us understand the origin of these singular- 
ites. They originate technically from the quark and antiquark propagators 
entering the Feynman amplitudes (the kernel ( |107|) is roughly their product): 

1 

(p + ky 
1 

{p' + kf 

The propagators diverge for two different configurations: 

i) CoUinear singularity, — > 0, i.e. when the gluon is emitted at a very 
small angle with respect to the source (the quark or the antiquark); 

ii) Infrared singularity. Eg 0, i.e. when the gluon is emitted with a very 
small energy. 

Let us see in detail how the end point singularites of (|107|) are related to the 
infrared and the collinear ones. The simple pole singularity 

Xi 1, Xj < 1 (117) 

implies 

Eg>0, 6,3 ^ 0, (118) 

i.e. the collinear emission of a gluon with a finite energy from the quark/ 
antiquark leg. The simple pole singularities ( |113| ) are therefore induced by 
collinears but not soft gluons. 
Let us consider now the infrared limit 

Eg ^ 0, e^gy^o, e^g^o (119) 

i.e. a gluon with small energy but not collinear to the quark or the antiquark. 
Kinematics simplifies for 

X3 < 1 (120) 



2EgEg{l - COS 6qg) S 1 " x 2 

1 1 1 



2EgEg{l — COSOgg) S 1—Xl 



(116) 
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to 

1 - COS 6qg 
Xi ~ 1 - X3 ^ 

1 - COS 9gg 
X2 — 1 — X3 — 

+ e-,g ^ n (121) 

such that 

COSOqg ~ — COS^^gg (l22) 

The infrared contribution to the total cross section is of the form 

f dxidx2 f dxidx2 
oa as 7 77 r ^ ots I ~t7 n — n — 7 (-'-23) 

J (1 — Xi)(l — 0:2) J Xi_[l — cos Oqg) [I — cos Oqg) 

Introducing the new variables 

e = X3 

u = cosOqg (124) 
so that ~ 1 — e(l + u)/2, X2 — 1— e(l — u)/2, we have in the new variables 

fa.^,/** (125) 

Note that the change of variables converts the illusory double pole in X3 into 
a simple one. Since u 7^ ±1, we have a simple pole singularity for e — > 0. 
Infrared singularities are therefore associated with simple poles. 

We see therefore that both coUinear and infrared singularities are associ- 
ated to simple poles of the differential cross section, i.e. to single logarithmic 
terms in the 3-jet cross section (see later). 

Let us consider finally the case (|112|) , i.e. when the energy fractions Xi 
reach simultaneously but independently their endpoint values. In this case 

Eg Q and 9qg ^0 or 9qg ^ Q (126) 

Eg. ( |125| ) is still valid, but in this case 

e ^ and m ^ ± 1, (127) 
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giving in this double pole singularity. We see that the leading singu- 

larity is generated by a gluon which is soft and at the same time coUinear to 
the source 

de rde 



5a ^ as I - j J (128) 



where we have taken m ~ 1 — ^^/2 in eg . (11251) . 



We can integrate the differential cross section ( p.07| ) as long as we avoid 
the endpoints. The related cross section will be finite but strongly dependent 
on the kinematical cuts. We note that these cuts are naturally introduced 
in an experiment in the form of resolution power. Since detectors have a 
finite angular resolution 5 > we cannot disinguish an isolated quark from 
a quark accompained by a gluon at an angle 9 < 6. Analogously, the finite 
energy resolution 77 > makes impossible to distinguish a quark from a quark 
surrounded by a gluon with energy e < 77. In other words, if the gluon is 
'too' collinear and/or soft, the final qqg is detected as a qq state: 

I qqg; e<r], e <5) \ qq)phys (129) 



The divergent total cross section (|115|) is therefore unmeasurable. The finite 



energy and angular resolution of the detectors restrict the integration region 
of ( |1071 ) for the observable qqg cross section, avoiding the singular regions. 
Technically: there is a non zero lower limit in the energy and angle integra- 
tions in ( |128| ). The observable cross section is therefore finite but dependent 
on the resolution parameters 6 and t]. 

In the same spirit, let us consider now the truly observable cross section for 

e+e" ^ qq (130) 

to order as- The qqg cross section (|107|) has to be integrated in the singular 
region, i.e. in the region corresponding to an undetectable gluon, but also 
virtual corrections to order 0:5 have to be included. It turns out by explicit 
calculation that collinear and infrared singularites cancel in the sum, giving 
rise to a finite cross section of order as- That is basically the same cancel- 
lation which occurs in the completely inclusive 0{as) correction (eqs.([77|) 
and (|79|)), which include the same virtual diagrams and the same singular 
integration regions for the real diagrams. 

To summarize, the divergences encountered mean that QCD has the ten- 
dency to produce gluons with small energy and with small opening angles 
with respect to massless sources. 
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CoUinear singularities (but not infrared ones) can be regulated with a 
small quark mass, which changes the propagators according to: 



{p + A;)2 - m2 2E,Eg{l - (3, cos O^g) 

(131) 



{p' + ky - m? '2EgEg{l - /3g COS 6- 



where (3 = p/E < 1 is the velocity. Therefore, the gluons emitted in the 
production of beauty quarks, i.e. in the process 

e+e" ^ b + b, (132) 

do not have any collinear singularity. The angular distribution of the gluons 
is given for small emission angles 6 by 

' (133) 



dcose l-pcosO l + (7^)2 

where 7 = E/m = — (3'^ ^ 1 is the Lorentz factor. The angular 

distribution has therefore a fast increase, 

dN ^ dN I , 
- (134) 

de dcose 9 ^ ' 

when we reduce the emission angle from ^ ~ 1 up to 

Orain ~ -• (135) 

7 

The behaviour ( |134| ) is the same as that of massless quarks. At even smaller 
angles 

e <emrn (136) 

the angular distribution is essentially constant. That is how the quark mass 
m acts as a collinear regulator: it cutoffs the increase of the angular dis- 
tribution. We expect therefore a logarithm of the mass after the angular 
integration: 

— ~ log— 137 

1/7 6 m 
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The angular distribution ( |133| ) is similar to the angular distribution of the 
classical electromagnetic radiation from a relativistic charged particle (see 

The above discussion about physical states and observable cross sections 
can be rephrased saying that we replaced single particle states with the true 
'm' and 'out' states, i.e. the states prepared with accelerators and observed 
with detectors. These are not anymore single particle states but have com- 
ponents with different number of particles 

I q)phys = a \ q) + b \ qg) + c \ qgg) + d \ qqq) + ... (138) 

The quark state is replaced by a quark surrounded by a cloud of collinear 
and/or soft gluons and massless qq pairs. In the next section we will give a 
simple prescription to build these truly observable out states. 

A remark is in order before ending this section. Up to now we considered 
the pure perturbative process of gluon radiation, so any angular or energy 
cut for the gluon is admissible (any y, see next section). The only limitation 
to the above results inside perturbative QCD is related to the Landau pole: 
as we lower the momentum transfer going into the soft region, the coupling 
rapidly increases leaving the perturbative region. But there is also a limita- 
tion of the applicability of perturbative QCD results related to confinement. 
A gluon with an energy of the order of the hadronic scale 

e ~ Aqcd (139) 

and/or with a transverse momentum with respect to the source of the 
order of the hadronic scale, 

kr ~ Aqcd (140) 

cannot move on a rectilinear motion neither producing a separate jet (not 
enough energy) nor flying away as an asymptotic state (confinement). This 
gluon is trapped by the color source. The perturbative behaviour of such 
gluon (an almost free motion) is therefore completely different from the 'real 
behaviour'. If we want perturbative QCD to describe the real hadronic 
world, we have to discard this conflict situation between perturbative and 
nonperturbative dynamics, imposing energy and angular cuts well above the 
hadronic scale: 

e > Aqcd, kr > Aqcd- (141) 
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That way we assign to different jets partons which can effectively produce 
different jets. 

5.2.1 The Jade Algorithm 

The definition of an N jet event according to the JADE algorithm is the 
following 0. Take the 4-momenta of the particles Pi and compute all the 
invariant masses: 

Mf, = ip^+PJ)^ 1^3 (142) 

If the smallest invariant mass (of particles i and j say) is less than a given 
threshold, 

M^. < ys, (143) 

then particles i and j are combined together to form a 'new' particle' (pseu- 
doparticle) with momentum pij = pi+pj. y <C 1 is a constant determining 
how 'fats' are the jets, such that ys ^ Aqcd to avoid large hadronization 
effects. If instead 

> ys, (144) 

particles i and j are not combined together. 

We iterate this procedure treating pseudoparticles as particles. At each 
step two particles are combined into a psedupoparticle, till the reduction 
stops because all the pair have invariant masses greater than the thresh- 
old. The final particles are identified with jets, and their number is the jet 
multiplicity. 

Experimentalists compute n jet fractions combining the meaured mo- 
menta of the hadrons in the final state, while theorists compuye n jet frac- 
tions with final states composed of partons. If the qualitative ideas about 
confinement and hadronization discussed in secs.^and ^]2|are correct, we can 
test perturbative QCD comparing hadronic jet fractions with partonic ones. 

In the case of a qqg final state, a 3-jet event is one in which all the three 
possible ivariant masses exceed the threshold: 

{p + kf > ys (145) 
{p' + kf > ys (146) 

{p + pT > ys (147) 
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which reads in terms of the energy fractions: 

x,<l-y (148) 

and therefore 

Xi + Xj>l + y ij^j (149) 

We see that the boundary in the integration domain (the dangerous terms) 
are ehminated. 

Let us define the jet fractions fi as 

J2 — 

O-tot 

h = (150) 

where atot = (^2jet + o'sje* = o'o(l + Q^s/^r) is the total cross section to order 
as (to compute to 0{as) we divide simply by ctq). It holds clearly: 

/2 + /3 = 1. (151) 

We have for /s 

„ _ /• , _ ^ "s [^-y dxi f^-y dx2{xl + xl] 

J3 — U.'^qqg — — 



27r J2j/ 1 — Xi Jl+y-xi 1 — X2 

( (^]h3- 6,) log ( -JL-] + 2 log^ ( » 



27i\ \l-yj °\l-2y) ° \^-y) 

where Li2{z) is the dilogarithmic function (also called the Spence function), 
defined by the integral representation: 

Lz2{z) = - r du (153) 
Jo 1 — u 



We see that there are no infrared singularities in eq. (|152D (like the 1/e poles 
we hit in eq.(|7^)): the jet fraction is an 'infrared safe' quantity. The jet 
fractions /s and /2 = 1 — /s are plotted in fig. 8 as a function of y and in fig. 
9 the experimental fi are represented. 
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If we vary the center of mass energy y/s of the e^e" coUision keeping 
y fixed, the jet fractions /« change because of the scale dependence of the 
couphng: 

fi{s,as,y) = fi{as{s),y). (154) 

We may therefore observe the variation of the QCD couphng measuring jet 
fractions at different energy with the same cut y. This is also a method for 
measuring as (see the lectures of B. Mele at this School). 



5.3 Sudakov Form Factors 

Let us assume that we have an experiment with a good resolution power in 
the measure of invariant masses, i.e. that we can meaure very small invariant 
masses. We can therefore measure the 2-jet fraction /2 = 1 — /s up to very 
small y 

y < 1 (155) 
The condition (|155|) is a strong restriction on the phase space of the final 



hadrons/partons: the final state consists of two very 'thin' jets in the space, 
which contain almost all the energy of the event. 

Neglecting terms which are infinitesimal with y in eq. (|152D (of the form 
y, 2/ log?/, ylog^y, etc., the so called power suppressed corrections), /2 reads 
to order as'. 

My) = l-CF^(2\og'y + 3hgy-^ + ^ + 0{y)j (156) 

There are three terms. There is a double logarithmic term which is the 
remnant of the collinear times infrared singularity, a single logarithmic term 
which is the remnant of the infrared + collinear singularity, and a finite term. 
Simbolically: 

fi^^ ~ aslog^y + aslogy + as (157) 

The y cut, as discussed previously, regulates both kind of singularites. If 
we use instead two different regulators for collinear and infrared singularities 
(for example an energy and an angular cut), we have for /2 an expression of 
the form: 

/i^^ ~ asLl + asL + asl + as (158) 
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where L is the colhnear logarithm (containing the angle cut) and / the infared 

logarithm (containing the energy resolutiom). 

The 0(q;|) corrections to /2 have also been computed [^: 

= {^)'{clZc{y) + CFNZM{y) + TnZT{y)^ (159) 

where: 

/13 \ /9 \ 

Zc{y) = 21og^|/ + 61og3 2/+(y-C(2))log2y+(--3C(2)-12C(3))logy 

+ ^-|c(2) + llC(3) + 4C(4) 

ZN{y) = ylog^i/+(2C(2)-^)log2y+(6C(3)-^)log|/ 

+ I + f C(2) - 13C(3) + ^C(4) 

ZAy) = -^log3i/ + ylog2y + 51ogy + ^-^C(2) (160) 

We see that the 0(a|) correction to /2 contains up to four logarithms of y. 
In general the correction to order ag contains terms up to log^" y, i.e. there 
at most two logarithms per loop (one colhnear times one infrared log for each 
loop). We have therefore: 

2n 
fc=0 

The perturbative expansion for f2{y) is not simply an expansion in powers 
of as at the relevant scale, like in the case of R. Every order in as is 
a polynomial in logy. That means that there are logarithms which are not 
absorbable in the renormalization of the coupling constant, as it does happen 
instead with R. Let us consider the physical difference between these two 
observable. R is an inclusive quantity, characterized by a single mass scale, 
the center of mass energy ^/s, so the theory can develop only logarithms of 
the form ^ 

log^, (162) 

s 

where Aq is the ultraviolet cutoff, which are absorbed by renormalization 
of the lagrangian since the theory is renormalizable. The jet fraction /2 is 
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instead a seminclusive quantity and its definition involves also another scale, 
ys. As a result, the theory can develop (and it actually does) logs of the 
other possible mass ratio, 

s 1 
log — = log - (163) 

ys y 

These logs are final, i.e. they do not cancel in the final result, as it does 
happen with the virtual and real diagram infrared singularities in R. The 
logs ( |163| ) have not ultraviolet origin (they are indeed infrared, as we saw), 
so they are not absorbable by a coupling redefiniton. 
For ?/ <^ 1 these logarithms can become so large 

log- > 1 (164) 

y 

that 

as\og^y ~ 1 (165) 

even though 

as < 1 (166) 

The convergence of the expansion in powers of as may therefore be spoiled by 
such large logarithmic coefficients. The solution of this problem is to abandon 
fixed order perturbation theory, in favour of an expansion according to the 
degree of singularity of the terms in the limit 

log- oo. (167) 

y 

The most singular terms, corresponding to the lowest order in this new 
asymptotic expansion, are of the form 

alloi^y n = 0,l,2,3,...,fc,.... (168) 

The sub-leading terms are of the form 

aglog^"-^?/ n = 1,2,3.. .fc... (169) 

and are of order l/logy <^ 1 smaller than the leading ones. The sub-sub- 
leading terms are of the form 

agW-^y n = 1,2,3.. .fc... (170) 
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and are of order 1 / log y smaller than the sub- leading ones and so on. Keep- 
ing only the leading series if often called double logarithmic approximation 
(DLA). Let us see how the two series are related. By simple rearrangement: 

= 1 + a^ax^ + aiax + + h^a^x'^ + hia^x^ + h2a^x^ + b^a^x + b^a'^ + . 
= (1 + a^ax"^ + boa'^x^ + ...j + (^aiax + bia^x^ + + {a2a + 62«^a;^ + ., 

+ {b^a''x + ..)+ ... (171) 

where x = logy, a more confortable notation for the coefficients has been 
used and we dropped the subscript on as. 

Since the coefficients of the leading terms ( |168| ) are known for any order 
(they exponentiate, see later), we can extend the validity of our result toward 
smaller y values factorizing the leading series. Up to order order we have 

/2 = (1 + aoax"^ + b^a^x'^ + ...^ {l + aiax + + (61 — aQai)a^x^ + 

(62 — aoa2)a^x^ + b^a^x + b^a'^ + ...) (172) 



We factorized the leading terms subtracting the spurious terms in the re- 
maining factor up to 0(q;^) included. The difference with the unfactorized 
expression is therefore 0{a^). In the first bracket on the right hand side we 
can add all the higher order terms of the leading series, i.e. resume the whole 
leading series, without affecting the second bracket up to order a^: 

/2 = (1 + a^ax^ + boQ^x"^ + cqo^x^ + ... + /oa^x^" + ... ) (1 + aiax + a2a 



+ (61 — aoOi)a^x'^ + (62 — ao«2)«^a;^ + fesa^x + 640^ + ... j (173) 

We may factorize the subleading series in the same way: 

/2 = {lead.)(j- + aiax + b[a'^x^ + ...^ (^l + a2a + b2a^x'^ + (63 — aia2)a^x + 64a;^^ 

(174) 

where b'l = bi — a^ai, b'2 = 62 — cto'^2 are the new coefficients shifted by the 
factorization of the leading series. This procedure can be carried on to any 
prescribed order in as- 

Note that the coefficient of the log^y term in ( |159| ) is 1/2 of that of the 
log^ y term in (|156|) : 
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/2 = l-C.^log^, + l(^) log^, + ... (175) 
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It can be proved that the leading terms of any order do exponentiate, i.e. 
they give a contribution to /2 of the form 



/2 = e-"sCFAiog y (^^7g) 

Expressions of this form are known as Sudakov form factors, who studied 
similar problems in QED long ago [110, ITR O, 123]. Since in lowest order 



/2 = 1, we see that the leading double logarithmic terms suppress the 2- 
jet cross section for y <C 1. There is a physical explanation for that. For 
y <^ 1 we are considering the quasielastic production of qq pairs at high 
energy, i.e. a final state with little activity around the tree level produced 
qq pair. Accelerated color charges naturally produce radiation, as it happens 
in classical electrodynamics. The Sudakov represents the suppression of the 
improbable non radiative channels. We see here an implementation of the 



factorization idea discussed in sec. 2.1: the basic electromagnetic process 



is corrected by soft gluon effects which appear as a factor in eq.( |17(j| ). The 
important physical information given by the Sudakov form factors is how 
the measured cross sections depend on the cut and the resolution of the 
experiments. 

Another important property of the Sudakov form factors is the broaden- 
ing of the sharp structures, like for example peaks in cross sections around 
resonances. Let us consider a QED case, the Z line shape discussed previ- 
ously. The Sudakov suppression of non radiative channels implies a relative 
increase of the radiative cross sections, so that energy is frequently released 
from the e+e~ system to the radiation field, through multiple photon emis- 
sions: 

e^e~ e~^e~ + n^y (177) 

The fluctuations in the energy released by the e"'"e~ pair induce equal fluc- 
tuations in the energy available for the production of the resonance, so the 
latter is not anymore produced always on the peak, for any selected beam 
energy. That smears the peak and shifts its position toward higher center of 
mass energies a/s > Mz- 

Let us observe that also the subleading terms proportional to [Cp)"', i.e. 
the QED-like ones, seem to exponentiate: 

h = l-%^(21og^, + 31og,) + if^)'(41og^, + 12W,)+... 



(178) 
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They do indeed exponentiate, producing the correction factor 

5/2 = e~^"^'"^/*-^'^''^°^^. (179) 

The terms (|1 79|) however do not contain the whole subleading correc- 
tions, as is clear looking to the form of the Z^ly) and Zriy) contributions 
in ( |160|) (they both contain a^log^y terms). The factorization of all the 
subleading corrections is a complicated problem because it involves many 
different phenomena: the variation of as with the scale (if for example 
a sin) asivfi) +0{allogy), 5{aslog^y) = 0(a|log^?/), i.e. a subleading 
term), separate coUinear and infrared singularities. 



6 Conclusions 

In these lectures we discussed the foundation of perturbative QCD and an im- 
portant physical application, the high-energy e~^e~ annihilation into hadrons. 
Even though perturbative QCD is technically similar to the perturbative ex- 
pansion of QED, the physical content is very different. In QED the physical 
one-particle states (electrons and photons) are identical to the ones in the 
free theory (free electrons and photons) after renormalization. One then 
makes a perturbative expansion of the scattering matrix elements, so only 
a (weak) convergence assumption of the series is postulated. In perturba- 
tive QCD the one-particle states (the hadrons) are not identical at all to the 
one-particle states in the free theory (free quarks and gluons), after renormal- 
ization. We neglect this problem and compute cross sections with quarks and 
gluons as they were asymptotic states. We relate subsequenty these unphys- 
ical processes to the physical ones with a set of qualitative ideas about their 
connection (parton-hadron duality). Apart from convergence problems, per- 
turbative QCD computations therefore have to be supplemented by assump- 
tions about non-perturbative phenomena (such as confinement, hadroniza- 
tion, jet formation), necessarily accompaining the perturbative process. In 
these lectures we discussed this point considerably, presenting also qualitative 
discussions and some models of confinement and jet formation. 

In the second part of the paper we considered the e^e~ annihilation to 
hadrons at high energy, which has many different dynamical aspects. The 
simplest observable is the total hadronic cross section. Its perturbative com- 
putation relies only on the assumption that hadronization does not change 
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the probability of the partonic process. Furthermore, the perturbative ex- 
pansion is very simple: it involves only numerical coefficients times power 
of as evaluated at the center of mass energy. More 'delicate' observables 
are the jet fractions. Their perturbative computation requires an extra as- 
sumption with respect to the total cross section. We have to assume local 
parton-hadron duality, i.e. that hadrons arc formed by partons which are 
close in phace space. This postulate is necessary because jet fractions con- 
tain dynamical effects related to intermediate scales, i.e. momenta p between 
a cutoff scale A ^qcd) and the center of mass energy 



These scales have not to be mixed by hadronization. Also the perturbative 
series for the jet fractions is more complicated: there are large logarithmic 
corrections of the form 



The latter is anyway a technical problem that is solved within perturbative 
QCD with resummation techniques. 

Comparing with experimental data, we have seen that the perturba- 
tive QCD approach, perturbative computations -|- hadronization assump- 
tions, does indeed work, both qualitatively and quantitatively. A face of 
the hadronic world seems to be largely understood. This also implies that 
confinement is very 'hidden' in many processes. The progress in the pertur- 
bative QCD direction consists in higher order computations, better resumma- 
tion techniques, and in finding measurable quantities that characterize some 
properties of the hadronic states and are little influenced by hadronization. 

A Evaluation of the 3-jet cross section 

In this section we compute the cross section at order as for 



< < s. 



(180) 





(182) 



We take the gluon propagator in Feynman gauge: 




(183) 
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and the quarks and leptons massless. The Feynman amphtude M is: 

M = Ma + Mb^ (184) 

where p — ^fi'if, s — q^, q — I + I' and I and I' are the 4- momenta of the 
electron and the positron. 

We take now the square of the modulus, average over the initial helicities 
and sum over the final helicities and polarizations, 

\M (185) 

rr'ss'h 

using the formulas 

^M^(p)m^(p) = p, T'^r{p')Vr{p') = P' 

r r 

Y.^H,{k)eUk) = -9,^ + ^ (186) 

(the k^j^ki,/ k"^ term in the photon polarization sum does not contribute because 

the diagrams are QED-like). 

We multiply by the relativistic final states factor 



d^p d^p' d^k 



(27r)32po (27r)32p'o {27r)V 
the relativistic normalization factor for the initial state 

1 1 1 



(187) 



;i88) 



the 4-momentum conserving 5 function 

{27iY6^''\q-p-p' ~k), (189) 

and divide by the relative velocity Vrei of the e^e~ pair in the center of mass 
frame (i.e. by the flux for volume = 1): 

Vrel = ^ = 2, (190) 
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(the two particles have opposite hght velocity). We have: 
da — 

Cpe'g^ 1 



^ 5: \M\''5^^\q-p-p'-k)^'^'^ 

Po P ^0 



f rr'ss'h 



(191) 



where the leptonic and the (completely integrated) hadronic tensors are given 

by 



d?p d?p' cPk 
Po P'o ko 



5^^\p + p' + k-q)H^,{p,p\k) (192) 



and i^ujy is the (unintegrated) hadronic tensor 



Tr 



J p + k p' + k \ p + k „ n p' + k 



\p + kf 



'{p' + kY 



(p + ky 



{p' + ky , 

(193) 

The tensors are symmetric under exchange of the indices \x and v and are 
transverse with respect to because of electromagnetic current conservation 



T^- = 0, q, L^^ = 0, 



Tfjtt, can be parametrized as 



so that 



The contraction gives 



pa- 



(194) 

(195) 
(196) 

(197) 



We decoupled the tensors with the projector guvQpal^- Therefore we need 
only to compute the contraction 



g^-" H,,{p,p',k) = 
Tr 



(198) 



q-p 



{q-p'f 



Q-P \-, a Q-P p p Q-P a 



{q-pf 



[q-p'f 



{q-pf 
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Because of the d^'^^ {p + p' + k — q) we made the replacements 
p + k ^ q—p\ p' + k ^ q — p. 

We have: 



+ 



{{q-p'W {(i-pf{(i-p'Y 



+ 



21 



+ 



22 



{q-pfiq-p^Y {{q-pY) 



where 



^11 


= +Tr 


P7p(9 


-p')ltiP'l^{q 


- P')Y 


-^^12 


= -Tr 




- p')inP'Y{q 


-p)r 


-^^21 


= -Tr 


P7m(? 


-phpp'i^iQ- 


-p')Y 


-^^22 


= +Tr 


P1h{q 


- phpP'YiQ - 


-p)Y 



The exphcit computation gives: 



X 



11 



4 Tr 



pqp' q 



X 



12 



-^^21 
-^^22 



4Tr[piq-p')p' {q-p') 
16 (^2p ■ q p' ■ q — q^ p ■ p'^ 

2 Tr P7p(g — ~ p)Yp' = ~ p) " ~ p') Tr 
32 p ■ p' {q — p) ■ {q — p') 

Xn 



P 



where we used the identities 



7/.7p7'' = -27p 



Tr 



IpHvlplu 



9p<T9yp 9ua9i^p 



ipncnaipi^ = -'^ipiisia- 
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Summing all the terms we derive: 

H^, = 16 (2p-qp' -q-q^p-p' 
p ■ p'{q -p)-{q- p') 



+ 



{{q-pfY {{q'-pY) 



+ 64 



{q-pYiq-p'f 



(204) 



The scalar projection of H^^ is much simpler to compute than the original 
tensor because it is easily reduced to the trace of four gamma matrices, while 
the computation of H^i, involves the trace of six gamma matrices. 

We perform now the integrations in the center of mass frame, where 



0. 



(205) 



First we integrate the 5^^-* {p-\- p + k) over the gluon 3-momentum which 
gives 

k ^ -p-0, (206) 

and therefore 



ko 



p + p 



PI + p'l + 2poP'o cos Oqq. 



(207) 



We choose now the spatial frame in such a way that p is directed along the 
z-axis and integrate over the polar angle of p'. 



P'o 



p'odp^d cos 9' d4>' , 



(208) 



The energy-conserving S{ko + po + Pq — go) -function gives: 
PodpodflpQdp'Qd4>'d cos 9' 



ko 



S{qo-po-Po-ko) g^^.H^'' ^ J dpodQdp',d<l>' g^^H^'' 

(209) 



because 



dkn 



1 -1 



d cos 9' 



kp 
PoPo 



The integration over dVLdcj)' is trivial and gives Stt^: 

dpodpodndcp'gf.^H"'' = 87r^Eldxidx2g^uH'"' 



(210) 



(211) 
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We express now the hadronic trace ( p04| ) in terms of energy fractions Xi and 
X2- We divide numerators and denominators by E^, and replace the rescaled 
scalar products according to the relations 



p-q = 2xi, p-q' = 2x2, 
These substitutions give: 



4, p-p' = 2(xi + X2-l). (212) 



1 — Xi 1 — X2 2{xi + X2 — 1; 

1 — X2 1 — Xi 



:i-xi)(i-x2) 



{1-Xi){l-X2) 



(213) 



We note that the leading singularity (1 — a;i)~^(l — 2:2)"^ in Feynman gauge 
comes from the interference term (while it comes from a single direct term 
in a particularly choosen axial gauge |Tl|, p!6[]). 

Putting all the pieces together, we have the final result: 



da 



2Cf a^as 



dxidx2 



^2 I /y»2 



(1 - Xi)(l - X2) 



(214) 



which coincides with result ( |107D multiplied by ctq = 47r/3a^/s. Performing 
the angular integrations and the Dirac algebra in n dimensions, we derive 
the dimensionally regularized cross section in eq.([77|). 
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FIGURE CAPTIONS 

Fig.l: A two jet event from hadronic Z decay (DELPHI, taken from ref. 
Fig.2: A three jet event from hadronic Z decay (DELPHI, taken from ref. 



Fig. 3: Born diagrams for e^e qq. 

Fig.4: Angular distribution of two jet events at 34 GeV center of mass energy. 
The curve is the QED + parton model prediction where x = cos 6 

(TASSO, from ref. |2|). 



Fig.5: Compilation of R values at low energy (upper figure) and at high energy 
(lower figure), (from ref. 



Fig.6: Total cross section for e^e^ +e.m. radiation, around the 

Z peak. The dotted line is the Born approximation, the dashed line 
is the Born approximation + leading logs coming from initial state 
radiation and the solid line is the Born approximation + leading logs 
+ sub leading logs (from ref. p6|). 

Fig. 7: Diagrams for real corrections (upper figure) and virtual corrections 
(lower figure) of order as to e'^e~ — >• qq. 

Fig.8: Plot of the jet fractions /2 and /s to order 0:5 (from ref. [p!2| ). 

Fig.9: QCD fits to the jet fractions (OPAL, from ref. pl). 
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